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Electron transport limited by the rotating exchange- 
potential of domain walls is calculated in the ballistic limit 
for the itinerant ferromagnets Fe, Co, and Ni. When realistic 
band structures are used, the domain wall magnetoresistance 
is enhanced by orders of magnitude compared to the results 
for previously studied two-band models. Increasing the pitch 
of a domain wall by confinement in a nano-structured point 
contact is predicted to give rise to a strongly enhanced mag- 
netoresistance. 
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The application potential of magnetoresistive effects 
has rekindled interest in the study of electrical trans- 
port in metallic (Stoner) ferromagnets such as Fe, Co, 
and Ni. One complicating factor which is still an open 
problem is the influence on the transport properties of 
the magnetic domain structure. Domain walls (DWs) 
result from minimizing the sum of the magnetostatic, 
magnetic anisotropy, and exchange energies and they 
can be driven out of a material by applying a mag- 
netic field. This modifies the transport properties but 
the magnitude and even the sign of the magnetoresis- 
tance MR = [i? sat - R ]/R = [G - G sat ]/G sat (where 
i?o = I/Gq and R sat = l/G sat are the zero- field and 
saturation-field resistances, respectively) remain a mat- 
ter of controversy. 

Early experimental studies on very pure iron samples 
showed a large MR of up to 90% at low temperatures 
jl| which was attributed to percolation through numer- 
ous domains B. Using a free electron model, Cabrera 
and Falicov [pf interpreted transport through a single 
DW as a tunneling process and the corresponding MR 
was found to be exponentially small. More recently, 
Tatara and Fukuyama Q calculated the DW conductiv- 
ity in the clean limit where the mean free path £ result- 
ing from defect scattering is much larger than the wall 
width Xdw- For a free electron model in a semiclassi- 
cal approximation they found an MR which scales like 
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—now/i^Dwkp), where udw is the domain wall den- 
sity and Ep = h 2 k F /2m is the Fermi energy. This MR 
is also very small for DW widths and Fermi energies of 
transition metals. In room temperature measurements 
of transport through Ni and Co films exhibiting stripe 
domain structures, Gregg et al. M measured significant 



negative MRs, much larger than predicted by any of the 
above theoretical work (but in a regime where £ < Xdw)- 
Levy and Zhang |j) subsequently pointed out that spin- 
dependent impurity scattering can strongly enhance the 
(negative) DW-MR. Breaking of the weak localization 
quantum correction by the exchange field leads to a pos- 
itive MR at low temperature Otani et al. and 
Riidiger et al. § measured a positive DW-MR in thin 
magnetic wires, but up to high temperatures. 

In the spirit of previous work on the giant magnetore- 
sistance of magnetic multilayers || we study DW scatter- 
ing in the ballistic limit, i.e. in the limit where the defect 
scattering mean free path £ is sufficiently larger than the 
system size {TCfl . These results are appropriate for clean 
point contacts with diameters d sufficiently smaller than 
£. We disregard lateral quantization, assuming d 3> Xf, 
with Xf the Fermi wavelength. For perpendicular trans- 
port in multilayers, the ballistic MR is of the same order 
of magnitude as the MR in diffuse systems |j 11 1 . When 
£ ^> Xdw, our calculated transport coefficients can be 
used as boundary conditions in the semiclassical Boltz- 
mann equations ]1^ | . 

The conductance G is given by Landauer's formula as: 



(i) 



where kn is the conserved Bloch vector parallel to the 
DW. The transmission amplitude of an incoming state 

\i to an outgoing state kit v through the DW sandwiched 
by the two bounding domains of the ferromagnet is de- 
noted by t v ^(k\\), k\\jjL and k\\v labeling flux- normalized 
states at the Fermi energy to the left and right of the 
DW, respectively, including the spin labels. 

A constant modulus for the local magnetization vector 
is assumed; its direction may be represented by a single 
rotation angle since we disregard the spin-orbit inter- 
action. 6 varies along the z-direction but is constant in 
the x, y-directions. The exchange field of the DW can be 
diagonalized by a local gauge transformation at the cost 
of an additional spin-rotation energy. Instead of treating 
this term by perturbation theory we employ here 
the WKB approximation, which has the important ad- 
vantage of being valid also for vanishing exchange split- 
tings. 

In order to understand the basic physics, let us 
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first consider a simple two-band model in which the 
plane wave states with parallel wave vector fc|| and en- 
ergy E^i — h 2 k 2 /2m ± A are modified by the do- 
main wall. In the WKB method the spinor wave 
functions are multiplied by z-dependent phase factors 
exp [if dz'J2mE ± (q(z'))/fi 2 - Jfejj]. The eigenenergies 
of the local Hamiltonian in which the gradient q(z) — 
88/ dz is taken to be constant are those of a "spin spi- 
ral" [pi: 



;2 r 



E±{q) 



2m 



k 2 + q 2 /A ± ^/k 2 q 2 +p 4 



(2) 



with A = h 2 p 2 /2m and k z determined by E±(q) 
I .- . The WKB-factor is imaginary for states propagat- 
ing through the whole DW and exponentially damped 
otherwise. In our adiabatic approximation we disre- 
gard all tunneling states, which is allowed in the limit 
AflH/fcf 3> 1. The eigenstates are not pure spin-states: 
the DW/spin spiral system acts like a spin-orbit scatterer 
to mix the two spin-directions. The DW conductance is 
thus limited by the local band structure with the small- 
est number of modes at the Fermi energy which is at the 
center of the DW where q is maximal, g m ax = k/^dw- 
For perpendicular transport |14|: 



G(q) 




y/4 



for q 2 < 2p 2 



p 2 +p 4 /q 2 for q 2 > 2p 2 



(3) 
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This equation holds when k F > q /A+p , i.e. when both 
spin bands are occupied. Note that transport parallel to 
the spin spiral is much less affected by the DW: Gii (q) = 
G(0) + O(<z 4 ). 

In bulk transition metals in which q 2 -C p 2 , the 
DW-MR becomes MR = (G(g max ) - G(0))/G(0) = 
{■k /2\owkF ) 2 for DW independent of the exchange split- 
ting. Using Eq. (J^), the Fermi wave vectors for one con- 
duction electron per atom, and the experimental width 
of the DW, we obtain the numbers in Table 1 for Ni, Co 
and Fe. The effect appears to be very small and likely 
to be swamped by other magnetoresistive effects such as 
the anisotropic or ordinary MR. The reason is clearly 
the smallness of the kinetic spin rotation energy as com- 
pared to the exchange splitting, i.e. q 2 <C p 2 ■ The DW 
only slightly deforms the Fermi spheres, resulting in a 
tiny magnetoconductance. In transition metals, however, 
many bands at the Fermi energy are much closer than the 
exchange splitting. When spin-up and spin-down states 
close to the Fermi energy are (nearly) degenerate, a DW 
which gives rise to a repulsive interaction between them 
may push the bands away from the Fermi energy and re- 
duce the conductance. Realistic band structures must be 
used in order to evaluate the importance of these split- 
tings. To this end we carried out first-principles calcula- 
tions of defect-free DWs in Ni(fcc), Co(fcc), and Fe(bcc) 



within the local spin density approximation (LSDA) to 
density functional theory (DFT). The open character of 
the leads can be captured by the embedding Green func- 
tion technique based on the linearized-augmented 
plane wave method (LAPW) and the muffin-tin shape 
approximation for the crystal potential. The transport 
coefficients and the conductance of samples with arbi- 
trary stacking of atomic monolayers with non-collincar 
spins can be calculated with this method. The technical 
details of the method are given in Refs. {[?-] ]. 

In the adiabatic limit the DW may be represented by 
a spin-spiral, which can be computed using conventional 
band structure techniques by the generalized Bloch the- 
orem based on a combined translation and spin-rotation 
operator jl5| . For narrow DWs a "linear" model is more 
accurate, in which we calculate the transmission (numer- 
ically) exactly for a magnetization which is rotated by a 
constant rate q ma x in a finite region of width Xdw be- 
tween single domain leads |p^| . 

The results are summarized in Table I for the two mod- 
els just considered for (i) experimental bulk DW widths 
and (ii) for DWs of monolayer width, both for a total spin 
rotation ir. Note the large difference between the first- 
principles calculations and the two-band model. Fig. 1 
displays the width-dependent DW conductance as a func- 
tion of the magnetization rotation angle per monolayer 
for Fe and Ni, respectively. We observe a linear depen- 
dence, MR ~ —(/max, in clear contradiction of the two- 
band model (Eq. (§)). 

We can understand these features using perturbation 
theory. The spin spiral can be represented by an interac- 
tion Hamiltonian which contains two operators H 
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and Hfr2 ~ q z , respectively [||. The energy band struc- 
ture of the bulk ferromagnet and thus the conductance is 
modified by this interaction. In non-degenerate pertur- 
bation theory the first-order term corresponding to 
vanishes. The second order term due to and first 

(2) 

order term due to _ff int both contribute to the order of 
g 2 , which explains the leading term in Eq. (||). How- 
ever, in the presence of degeneracies simple perturbation 
theory breaks down. Instead, the Hamiltonian must be 
diagonalizcd first in the subspace of (nearly) degenerate 
states. The splitting of the degenerate states is directly 
proportional to the matrix elements of the interaction 
Hamiltonian, thus in leading order proportional to q. As 
the energy splittings increase, conducting channels are 
removed from the Fermi surface and the conductance is 
reduced proportionally. The linear dependence observed 
in Fig. 2 can thus be explained by the occurrence of many 
(nearly) spin-degenerate states close to the Fermi energy. 
Naturally, the MR is also much larger for closely spaced 
states which are not strictly degenerate. This explains 
the large difference between the results for the two band 
model and the full band structures in Table I. 

We observe that the relative effect of the DW is still 
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rather small, smaller than the experiments by Gregg jj| 
and smaller than the theoretical results by Levy and 
Zhang for very spin-asymmetric bulk defect scatter- 
ing p. Bloch DWs in thin films can be significantly 
narrower than in bulk material Jl9| , which means that 
the bulk DW magnetoresistance should be larger in thin 
films than the bulk values in Table I, but still smaller 
than in [^|J^]. The present calculations show unambigu- 
ously that the DWs increase the resistance. The experi- 
mentally observed DW-induced decrease of the resistance 
can therefore not be an intrinsic effect, but must be 
an as yet unidentified defect-related, size-related, or other 
extrinsic phenomenon. Previous results obtained by per- 
turbation theory and a two-band model |||| should be 
reconsidered in the light of the present findings. Unfor- 
tunately, implementing degenerate perturbation theory 
for diffuse systems 0] with realistic band structures ap- 
pears to be quite cumbersome E(J . In the recent work of 
Levy and Zhang Q the DW scattering is calculated on 
the basis of a two-band model. In spite of the small for- 
ward spin-flip scattering in this model they can explain 
a significant MR due to a strongly spin-dependent bulk 
defect scattering. The band structure crowding at the 
Fermi surface enhances not only the backward scattering 
which causes the ballistic MR discussed here, but also 
the forward scattering. The bulk defects might therefore 
be less important than initially apparent. 

The DW scattering increases with q = 89 /dz, which 
can be achieved by reducing the DW width or by increas- 
ing the winding number for a given width. Both oper- 
ations become possible by trapping a DW in a nanos- 
tructured ferromagnetic point contact. Ballistic point 
contacts have been fabricated successfully in simple met- 
als [^1 22 1, but not yet in ferromagnetic materials p3[ . 
When the magnetization on one side of the contact is 
pinned by shape anisotropy or exchange biasing, the mag- 
netization vector on the other side can be rotated inde- 
pendently by rotating the sample in an external mag- 
netic field. The maximal effect is expected for an abrupt 
domain wall, for which we predict a huge MR (see Ta- 
ble I) , much larger than what has been be achieved with 
tunnel junctions of the same materials. The material de- 
pendence on the angle between the two magnetizations 
(Fig. 2) betrays again the importance of the details of 
the band structure. In a similar fashion an tw-DW could 
be created by repeated rotation in the magnetic field. 
The conductance is then predicted to decrease linearly 
with the number of turns as in Fig. 1, up to some value 
at which phase-slips occur, or the spiraling magnetiza- 
tion spills out of the constriction. We stress that this 
somewhat naive picture needs to be supported by micro- 
magnetic calculations [ p3[ . 

In conclusion, we presented and analyzed model and 
first-principles calculations of electron transport through 
magnetic domain walls. The large number of bands close 
to the Fermi surface causes a strong enhancement of the 



DW-MR as compared to two-band calculations. Evi- 
dence that degeneracies at the Fermi surface of Fe, Co, 
and Ni can give rise to relatively large effects is found. 
DWs always decrease the ballistic conductance, causing 
a negative MR. The ballistic DW magnetoresistance is 
found to be somewhat smaller than measured recently 
in thin films, which can be partly due to the reduction 
of domain wall widths in thin films as compared to bulk 
ferromagnets. Trapping a domain wall in nanostructured 
constrictions is predicted to give rise to a strongly en- 
hanced magnetoresistance. 
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Ni 


Co 


Crystal structure 


bec 


fee 


fee 


Layer direction 


(100) 


(100) 


(111) 


Gsat [lO 15 !^- 1 ™- 2 ] 


1.531 


1.923 


1.529 


DW thickness Xdw 








(in nm) 


40 


100 


15 


(in monolayers) 


276 


570 


72 


Spiral angle/monolayer 


0.65° 


0.32° 


2.5° 


DW-MR 








2-band model 


-0.0008 % 


-0.0001 % 


-0.008 % 


adiabatic model 


-0.13 % 


-0.11% 


-0.33% 


linear model 


-0.39 % 


-0.16% 


-0.46% 


abrupt DW 


-71 % 


-58 % 


-67 % 



TABLE I. Parameters for Fe, Ni and Co, calculated satu- 
ration (single-domain) conductances and magnetoresistances 
(MR) as defined in the text. DW thicknesses are taken from 
Ref. §J. 
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FIG. 1. Conductance of domain walls in Ni and Fe as a 
function of the magnetization rotation angle per monolayer, 
A9 = hcl/\dw , where a is the monolayer width and \dw the 
width of the domain wall. Results are given for the adiabatic 
approximation (spin spiral) and the linear approximation (see 
text). The bulk ballistic conductances are indicated by the 
horizontal lines. 



FIG. 2. Conductances of abrupt domain walls in Ni, Fe, 
and Co as a function of the angle A(j> between the magneti- 
zation vectors of the bounding domains. 



than a constant chemical potential approximation. It can 
be shown that the correction to the conductance due 
to the shift of the chemical potential is SG(q)/G(0) = 
— g 2 p 4 /48fc F to the lowest order in A/Ef, which is neg- 
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